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1 . INTRODUCTION 


The  computer  code  EMPFIT  provides  a method  for  fitting  an 
electromagnetic  pulse  (EMP) , as  well  as  other  traces,  with  a simple 
function  that  is  easily  differentiated  and  Fourier  transformed.  This 
report  is  a supplement  to  the  original  report  on  EMPFIT  and  reflects 
the  changeover  from  a Control  Data  Corp.  (CDC)  6600  computer  at  the 
Mobility  Equipment  Research  and  Development  Command  (MERADCOM) , Fort 
Belvoir,  VA,  to  the  IBM  370/168  computer  at  the  Harry  Diamond 
Laboratories  (HDL) , Adelphi,  MD.  Along  with  this  conversion  effort, 
several  modifications  have  been  added,  making  EMPFIT  more  useful  and 
versatile. 


In  the  following  sections,  a description  of  the  time  derivative 
calculations  and  options  detailing  read  and  write  selections  on  a disk 
file  of  EMPFIT  will  be  given.  Next,  a section  on  increasing  the 
accuracy  of  the  code  using  the  AUTODBL  option  of  the  IBM  370/168  is 
explained.  Also,  a time  cutoff  modification  will  be  detailed  and  a full 
listing  of  all  of  the  output  options  is  included.  A brief  description 
of  the  new  subroutines  added  to  EMPFIT  is  discussed  and  an  up-to-date 
section  on  preparing  input  cards  to  run  the  code  is  also  presented. 
Finally,  a section  on  the  job  control  language  (JCL)  necessary  to  run 
EMPFIT  at  the  Woodbridge  Research  Facility  (WRF)  is  given.  In 
appendices  A and  B,  a listing  and  sample  run  of  EMPFIT  along  with 
examples  of  input  data  appear. 


2.  DERIVATIVE  OF  FITTING  FUNCTION 


Following  the  theory  of  EMPFIT 
1 < i < N using  the  function 


lThomas  V.  Noon,  User's  Manual  for  the  Modular  Analysis-Package 
Libraries  ANAPAC  and  TRANL,  Harry  Diamond  Laboratories  TR-1782-S 
(September  1977) . 


3.  READ  AND  WRITE  DISK  FILE  OPTIONS 

Since  there  was  an  apparent  need  to  be  able  to  use  data  not  in  card 
format,  the  necessary  modifications  were  made  to  utilize  data  that  exist 
on  a permanent  disk  file.  With  these  new  changes,  input  data  to  EMPFIT 
can  now  be  entered  through  a disk  file  and  it  is  also  possible  to  add  or 
delete  specified  data  pairs  from  this  disk  file  data.  When  the  "read 
from  a disk  file"  option  is  chosen,  a numbered  printout  of  the  data 
from  the  disk  file  is  produced  (see  fig.  3)  and  any  deletions  or 
additions  are  made  in  reference  to  this  listing.  Deletions  are 
accomplished  by  selecting  the  integer  number  of  the  time -amplitude  point 
to  be  scratched  and  entering  this  value  as  described  in  section  8.  An 
example  of  the  input  cards  necessary  to  delete  data  points  appears  in 
appendix  B.  Additions  to  the  disk  file  data  are  entered  by  the 
time-amplitude  value  and  are  automatically  placed  in  the  proper 
location.  The  use  of  this  option  is  also  detailed  in  section  8,  while 
an  example  of  the  input  cards  required  to  add  data  points  can  be  found 
in  appendix  B. 


Figure  3.  Example  of  disk  file  printout 
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figure  3.  Example  of  disk  file  printout  (cont'dl. 
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Figure  3.  Example  of  disk  file  printout  (cont'd) 
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In  addition  to  being  able  to  read  input  data  from  a permanent  file, 
it  is  now  possible  to  write  specific  output  to  a permanent  disk  file  so 
that  other  computer  codes  may  use  EMPFIT's  calculations.  Thus,  EMPFIT 
is  now  set  up  to  write  the  calculations  of  the  derivative,  f*  (t) , on  a 
disk  file  at  either  the  input  data  points  or  the  curve  fit  data  pairs. 
Also,  the  calculations  for  the  curve  fit,  f(t),  may  be  written  on  a 
permanent  file.  The  utilization  of  these  permanent  file  disk  writes, 
which  are  performed  on  separate  devices,  is  fully  described  in  section  8. 


It  should  be 
data  from  or  to 
READPF  and  WRITPF 
write  operations, 
need.  Figures  4 
starred  cards 

particular  user’s 
requires  the  user* 
//FORT . SYS IN  DD  * 


noted  that  it  is  possible  in  EMPFIT  to  read  and  write 
a disk  file  in  any  particular  format.  The  subroutines 
are  used  by  EMPFIT  to  perform  the  disk  file  read  and 
respectively,  and  may  be  changed  to  fit  the  user's 
and  5 show  the  listings  of  READPF  and  WRITPF  and  the 
reflect  cards  that  may  be  changed  according  to  a 
requirements.  To  employ  these  user  modifications 
s version  of  READPF  and  WRITPF  to  be  placed  behind  the 
card  seen  in  section  9. 


SUBROUTINE  RE ADPF (NT  ,NSETS  ,X ,Y , JPTS J 
DIMENSION  X(JPTS1 , Y(JPTSI 
DO  10  J*1,NSETS 

READ  (NT)  N1,N2,JPTS,(XU),A,I*1,JPTS)  • ( Y ( I)  ,B. 1*1, JPTS) 
10  CONTINUE 
RETURN 
END 


Figure  4.  Listing  of  subroutine  READPF. 


SUBROUTINE  WRITPF (NT  ,X  ,Y  ,NPTS  * AT  IT LE ) 
DIMENSION  X(NPTS), Y(NPTS ), ATI TLE  < 10) 
WRITE(NT)  ATITLE  * 

WR)TE(NT)  NPTS  * 

WRITE(NT)  X t Y * 

RETURN 

END 


Figure  5.  Listing  of  subroutine  WRITPF. 


4.  DOUBLE  PRECISION  AND  EXPONENTIAL  CALCULATIONS 


The  IBM  version  of  EMPFIT  has  been  written  in  such  a manner  that  the 
AUTODBL  option  of  the  HDL  IBM  370/168  can  be  used  and  not  affect  the 
plotting  software.  This  variation  is  seen  in  the  subroutine  ANOTAT 
which  has  been  constructed  differently  from  the  CDC  version  to 
accommodate  the  AUTODBL  option,  which  automatically  converts  all  single 
and  double  precision  variables  to  double  and  extended  precision 
quantities,  respectively.  This  conversion  allows  greater  accuracy  in 
the  computations  and  is  easy  to  employ.  Only  minor  changes  are  required 
in  EMPFIT,  specifically  dealing  with  the  plotting  software  of  ANAPAC.* 
In  particular,  the  only  modifications  occur  in  calling  the  double 
precision  versions  of  subroutine  DRAW4 , which  appear  in  EMPFIT' s 
plotting  routine  PLOTT.  If  the  AUTODBL  option  is  chosen,  EMPFIT  should 
be  compiled  with  the  following  changes  in  subroutine  PLOTT: 


(1)  The  DRAW4  calls 

CALL  DRAW4 (1,3,3,3,8,20,  XTITLE , YTITLE,  ATITLE,  TITLE) 
CALL  DRAW4 (2,3,1 LNLOG,  IPTS , -2 , 10 , T, F , 0 . , 0 . ) 

CALL  DRAW4( 2,3, ILNLOG,  MAXPTS,  0 , 10 ,TT , FF ,0 . , 0 . ) 

CALL  DRAW4 ( 3 , 3 , 0 , 0 , 0 , MAXPTS,  TT,FF,2.,0.) 

should  be  changed  to 

CALL  DRAW41 (3,3,3,8,20,  XTITLE,  YTITLE,  ATITLE,  TITLE) 
CALL  DRAW42 ( 3 , ILNLOG , IPTS,  -2 , 10 ,T ,F , 0 . , 0 . ) 

CALL  DRAW42 (3 , ILNLOG,  MAXPTS,  0 , 10 , TT , FF , 0 . , 0 . ) 

CALL  DRAW43 (3, 0,0,0,  MAXPTS,  TT,FF,2.,0.) 

(2)  The  DRAW4  calls 

CALL  DRAW4 (1,3,3,3,8,20,  XTITLE,  YTITLE,  ATITLE,  TITLE) 
CALL  DRAW4 (2,3,  ILNLOG,  MAXPTS,  0 , 10 ,TT ,FF , 0 . ,0. ) 

CALL  DRAW4 ( 3 , 3 , 0 , 0 , 0 , MAXPTS,  TT,FF,2.,0.) 

should  be  changed  to 

CALL  DRAW41 (3,3,3,8,20,  XTITLE,  YTITLE,  ATITLE,  TITLE) 
CALL  DRAW42(3,  ILNLOG,  MAXPTS,  0 , 10 ,TT , FF ,0 . , 0 . ) 

CALL  DRAW43 (3, 0,0,0,  MAXPTS,  TT,FF,2.,0.) 


1 Thomas  V.  Noon,  User's  Manual  for  the  Modular  Analysis-Package 
Libraries  ANAPAC  and  TRANL,  Harry  Diamond  Laboratories  TR-1782-S 
(September  1977). 
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(3)  The  DRAW4  calls 

CALL  DRAW4 (1,3,4,5,14,20,  XTITLE,  YTITLE,  FTITLE,  TITLE) 

CALL  DRAW4 (2,3,2, OPTS , 0,10, OMEGA,  ZABS,  0.,0.) 

CALL  DRAW4 (3 , 3 , 0 , 0 , 0 , OPTS , OMEGA , ZABS , 2 . ,0. ) 

should  be  changed  to 

CALL  DRAW4 1(3,4,5,14,20, XTITLE , YTITLE , FTITLE , TITLE) 

CALL  DRAW42 (3 , 2 ,OPTS ,0 , 10 , OMEGA, ZABS , 0 . ,0.) 

CALL  DRAW43 ( 3 , 0 , 0 , 0 ,OPTS , OMEGA , ZABS , 2 . , 0 . ) 

(4)  The  DRAW4  calls 

CALL  DRAW4 (1,3,3,3,12,20, XTITLE , YTITLE , FTITLE , TITLE ) 

CALL  DRAW4 (2,3, ILNLOG ,MAXPTS ,0,10, TT , DFF , 0 . ,0. ) 

CALL  DRAW4 ( 3 , 3 , 0 , 0 , 0 , MAXPTS , TT , DFF , 2 . ,0.) 

should  be  changed  to 

CALL  DRAW4 1(3,3,3,12,20, XTITLE , YTITLE , FTITLE , TITLE) 

CALL  DRAW42 (3 , ILNLOG , MAXPTS, 0 , 10 ,TT , DFF , 0 . ,0.) 

CALL  DRAW43 (3, 0,0,0, MAXPTS, TT, DFF, 2. ,0.) 

It  should  also  be  noted  that  a convention  involving  the  exponential 
function  calculations  has  been  employed.  This  is  due  to  the  word  size 
difference  between  the  CDC  6600  computer,  which  has  60  bits  per  word, 
and  the  IBM  370/168  computer,  which  has  only  32  bits  per  word.  Since  the 
CDC  version  of  EMPFIT  allows  exponential  arguments  to  be  taken  up  to 
approximately  720,  while  the  IBM  exponential  argument  size  is  limited  to 
approximately  174,  a simple  change  was  added.  This,  change  was  mainly 
influenced  by  the  fact  that  EMPFIT  fits  data  primarily  in  the  time  frame 
of  1 to  5000  shakes  (1  shake  ■ 10_8  s)  and  exponential  values  sometimes 
need  to  be  calculated  which  exceed  the  IBM  limit  of  174.  Thus,  it  was 
decided  that  the  easiest  solution  was  to  limit  the  size  of  exponential 
values  to  174  rather  than  to  make  extensive  software  changes  to  EMPFIT. 
This  limitation  is  only  used  when  fitting  the  exponential  functions  (eq 
1 and  3)  to  the  front  and  rear  of  the  fitting  function  where  the 
values  for  a and  8 in  equations  1 and  3 are  minimized  so  that  the 
argument  size  does  not  exceed  174.  When  this  convention  is 
automatically  employed,  EMPFIT  prints  out  the  largest  values  that  can  be 
used  for  a.  and  0 and  the  message  that  this  variation  has  been  utilized. 


5.  TIME  CUTOFF  OPTION 

Since  it  is  sometimes  necessary  to  deal  with  data  with  amplitudes 
very  nearly  zero  (with  respect  to  the  peak  amplitude)  at  the  tail  of  a 
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waveform,  an  addition  has  been  employed  which  allows  the  user  to  chop 
off  the  trace  after  a prescribed  time  value.  This  option  allows  one  to 
ignore  meaningless  information  at  the  end  of  a waveform  which  exists  on 
a permanent  disk  file  and  cannot  be  deleted  easily.  An  example  of  this 
option  can  be  seen  in  figures  6 and  7.  The  use  of  this  modification  is 
detailed  in  section  8. 
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Figure  6.  Example  of  curve  with  long,  nearly  zero  tail. 
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Figure  7.  Example  of  time  cutoff  option  applied  to  waveform 
with  long,  nearly  zero  tail. 

| 

6.  OUTPUT  OPTIONS 

With  the  aforementioned  modifications  EMPFIT  now  consists  of 
numerous  output  options.  Thus  a total  compilation  of  these  selections 
will  now  be  listed.  The  use  of  these  options  is  fully  detailed  in 
section  8. 

a.  Multiple  runs 

b.  Plot  titles  and  axes  information  supplied  in  program 


(1) 

T 

VS 

(2) 

T 

vs 

(3) 

T 

vs 

(4) 

T 

vs 

(5) 

T 

vs 

(6) 

T 

vs 

Er — time  versus  electric  field  strength — radial 
Ev~ time  versus  electric  field  strength — vertical 
B^— time  versus  magnetic  induction — azimuthal 
JR — time  versus  current  density — radial 
Jv— time  versus  current  density — vertical 
o— time  versus  air  conductivity 
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Plot  titles  and  axes  information — may  be  entered  by  the  user  if 
information  supplied  by  the  program  is  not  sufficient 


Enter  input  data 

(1)  Read  data  points  from  cards 

(2)  Read  data  points  from  a disk  file 

(a)  Delete  points  from  the  disk  file  data 

(b)  Add  points  to  the  disk  file  data 

Number  of  points  calculated  between  input  data  points 
Maximum  time  to  be  calculated 

Cut  off  time  and  amplitude  values  after  a prescribed  time 
Calculate  the  time  derivative,  f ' (t) 

(1)  f ' (t)  calculated  at  input  data  points  only 

(2)  f’ (t)  calculated  at  curve  fit  points 

Calculate  the  Fourier  transform 

(1)  Starting  frequency 

(2)  Maximum  frequency  to  be  calculated 

Plots 

(1)  Linear 

(2)  Log- log 

(3)  Semilog  Available  on  curve  fit  and 

derivative  plots  only; 

(a)  Abscissa  (log)  Fourier  transform  plots 

(b)  Ordinate  (linear)  are  always  log- log 

(4)  Minimum  value  of  ordinate  to  be  plotted  (does  not  plot 
data  points  below  this  value) — useful  when  plotting 
linear  or  semilog  and  the  range  of  values  is  very  large 

(5)  Curve  fit  plots— three  plots  may  be  obtained: 


(a)  Data  points  (plotted  with  X's)  and  curve  through 
data  points 

(b)  Curve  only 

(c)  Both  (a)  and  (b) 
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(6)  Derivative  plot 


(7)  Fourier  transform  plot 


Write  data  on  a disk  file 


(1)  Write  the  curve  fit,  f(t),  on  a disk  file 

(2)  Write  the  time  derivative,  f ' (t) , on  a disk  file 


7.  DESCRIPTION  OF  NEW  SUBROUTINES  OF  EMPFIT 

The  modified  version  of  EMPFIT  contains  the  following  new 
subroutines. 

PERFIL  Subroutine:  reads  the  data  from  the  disk  file  and  writes 
the  data  points  on  output;  deletes  or  adds  any  specified 
data  points  to  the  disk  file  data. 

CUTOFF  Subroutine:  cuts  off  the  time  and  amplitude  values  after 
a prescribed  time  value. 


DERIV  Subroutine : calculates  the  time  derivative  at  either  the 
input  points  or  the  curve  fit  points. 

RITDAT  Subroutine:  writes  either  or  both  the  time  derivative  data 
and  the  curve  fit  data  on  a disk  file. 


READPF  Subroutine:  reads  input  data  from  a disk  file. 

WRITPF  Subroutine : performs  the  actual  writing  of  the  data  onto 
the  disk  file  as  specified  by  RITDAT. 

8.  DATA  INPUT  PREPARATION  FOR  EMPFIT  (REVISION) 

Input  data  cards  for  EMPFIT  are  prepared  in  the  following  manner. 
Examples  of  input  card  decks  appear  in  appendix  B. 


Column 


Variable 


Format  Explanation 


Card  1 : Multiple  Run  Card 


Number  of  runs 


17 


Column  Variable  Format  Explanation 

Card  2 < Plot  Parameter  Card 


10 

IDENT 

11 

Identifies  data  to  be  read  in 
IDENT  =1  T vs  E 

=2  Tvs  Ej 

=3  T vs 

=4  T vs  J 

“5  T vs  Jy 

“6  T vs  cr 

20 

IFFT 

11 

Fourier  transform? 

IFFT  =0  Yes 
*=  1 No 

30 

I PLOT 

11 

Plots? 

I PLOT  = 0 Yes 
* 1 No 

40 

ILNLOG 

11 

Plots  in  linear,  log-log, 
or  semilog? 

ILNLOG  « 1 Linear 

= 2 Log-log 

■ 3 Semi log 

Abscissa  (log) 
Ordinate  (linear) 

Note : Fourier  transform  plots  are  always 
log-log. 

41-50  ORDMIN  E10.3  Minimum  value  of  ordinate  to 

be  plotted  (all  points  below  this 
value  are  not  plotted  on  linear 
or  semilog  plots) 


60 

ICURV 

11 

Curve  fit 

plots 

ICURV  - 0 

One  plot  of  data 
points  and  curve  fit 

= 1 

One  plot  of  curve 
fit  only 

18 

= 2 

Two  plots— one  plot 
of  data  points  and 
curve  fit,  one  plot  of 
curve  fit  only 

t 

r 

r 

■ 1 

J 

A • 

[ i 

Column 

Variable 

Format 

Explanation 

[ 

70 

IOT 

11 

Enter  own  plot  and  axes  labels? 

t 

IOT  - 0 Use  plot  and  axes 
labels  supplied  in 
program 

■ 

1 

- 1 Enter  own  labels  on 
cards  6,  7,8 

, * 

80 

IREAD 

11 

Read  data  points  from  permanent 
file? 

IREAD  - 0 No 
- 1 Yes 

Calculate  time  derivative? 

h 

IREAD  - 0 No 
- 2 Yes 

.as 

Do  both  of  above? 

IREAD  - 0 No  » 

- 3 Yes 

St 

Card  3i  Title  Card 

: i 

I ! 

1-80 

TITLE 

8A10 

Title  or  subtitle 

Card  4 r Fitting 

Parameter  Card 

•i 

8-10 

IPTS 

13 

Number  of  data  points 

read  in  j 

18-20 

MPTS 

13 

Number  of  points  calcu- 
lated between  input 
data  points 

21-30 

TMAX 

E10.3 

Maximum  time  to  be 
calculated  in  curve 
fit  calculations 

31-40 

ALPHA 

E10.3 

Used  to  fit  f(t)  ■ 

Ajeot  + A3e2at  to 

front  of  waveform 
(good  starting  value* 
a - 1.2E+8) 

19 


t 


1 


Column 

Variable 

Format 

Explanation 

41-50 

BETA 

E10.3 

Used  to  fit  f(t)  - 
A2e“®t  + AL(e“2®t  to 
end  of  waveform  (good 
starting  value:  6 ■ 
5.0E+4) 

51-60 

OSTART 

E10.3 

Frequency  to  start 
Fourier  transform  calcu 
lations 

61-70 

OMAX 

E10.3 

Maximum  frequency  to  be 
calculated 

80 

ICUT 

11 

Time  cutoff? 

I CUT  = 0 No  time  cutoff 
= 1 Cutoff  time 

points  at  pre- 
scribed time 
shown  on  card  5 

Note:  If  IREAD  = 0 skip  to  card  13 
* 1 Skip  to  card  10 
= 2 Go  to  card  9 
= 3 Go  to  card  9 

If  IOT  = 1 Go  to  cards  6,7 
8;  then  go  to 
IREAD  prescribed 
cards 

If  ICUT  = 1 Go  to  card  5; 
then  go  to 
IREAD  and  IOT  pre- 
scribed cards 

Card  5:  Time  Cutoff  Card 

1-10  TCUT  E10.3  Prescribed  cutoff  time 

value 

Card  6:  Abscissa  Label  Card 

1-10  XTITLE  A10  X label;  start  in  column  1 


r* 


Column 

Variable 

Format 

Explanation 

Card  7s 

Ordinate  Label  Card 

1-20 

YTITLE 

2A10 

Y label;  start  in  column  1; 
on  output,  ordinate  label  is 
only  in  A10,  A2  format  instead 
of  2A10 

Card  8s 

Plot 

Label  Card 

1-40 

ATITLE 

4A10 

Plot  label;  start  in  column  1 

Card  9s 

Time 

Derivative  Card 

j 

10 

IFPRIM 

11 

IFPRIM  * 1 f ' (t)  calculated  at 
input  data  points 
only 

=2  f ' (t)  calculated  at 
input  and  calculated 
points 

20 

IPERM 

11 

Write  f ' (t)  on  permanent  file? 

IPERM  = 0 No 

= 1 Yes  (as  specified  by 
IFPRIM) 

Write  f(t)  on  permanent  file? 

I PERM  = 0 No 
= 2 Yes 

Do  both  of  above? 

IPERM  = 0 No 


Column 


Variable 


Format 


Explanation 


Card  10:  Permanent  File  Card 


8-10  NSETS  13 


8-20  IDELET  13 


Number  of  particular  data 
set  to  be  read  from  disk 
file 

Delete  points  from  disk  file 
data? 

IDELET  = 0 No 
= 1 Yes 


28-30  I ADD 


13  Add  points  to  disk  file 

data? 


I ADD  = 0 No 
= 1 Yes 

Note:  If  IDELET  = 1 Go  to  card  11 

If  I ADD  = 1 Go  to  card  12 
If  IDELET  & I ADD  = 1 Go  to  card  11, 

then  12 

If  IDELET  & I ADD  - 0 Input  for  this 

data  run  is 
complete 


Card  11 : Deletion  Card 


8-10 

NUM 

13 

Number  of  data  point  on 

18-20 

NUM 

13 

disk  file  that  is  to  be 

28-30 

NUM 

13 

deleted  from  disk  file  data 

38-40 

NUM 

13 

48-50 

NUM 

13 

58-60 

NUM 

13 

68-70 

NUM 

13 

78-80 

NUM 

13 

Note : At  present  only  eight  data 
pairs  can  be  deleted. 


Column 

Variable 

Format 

Explanation 

Card  12: 

Addition  Card 

8-10 

NIADD 

13 

Number  of  data  pairs  to  be 
added  to  disk  file  data 

11-20 

TX 

E10.3 

Time  value  of  point  1 to  be 
added 

21-30 

FY 

E10.3 

Amplitude  value  of  point  1 
to  be  added 

31-40 

TX 

E10.3 

Time  value  of  point  2 to  be 
added 

41-50 

FY 

E10.3 

Amplitude  value  of  point  2 
to  be  added 

51-60 

TX 

E10.3 

Time  value  of  point  3 to  be 
added 

61-70 

FY 

E10.3 

Amplitude  value  of  point  3 
to  be  added 

Note:  Additional  cards  may  be  entered 
starting  in  column  11.  Data  input  is 
complete. 

Card  13: 

Data 

Card 

1-10 

T 

E10.3 

Time  value  of  first  point 

11-20 

F 

E10.3 

Amplitude  of  first  point 

Card  14: 

Data 

Card 

1-10 

T 

E10.3 

Time  value  of  second  point 

11-20 

F 

E10.3 

Amplitude  of  second  point 

Note:  Card  13,  card  14,  . . . are 
repeated  with  respect  to  the  number 
of  input  points  indicated  on  card  4, 
IPTS. 

Note:  Card  2 to  card  13,  card  14, 

. . . are  repeated  according  to  the 
number  of  times  identified  on  card  1, 
NRUN. 
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9.  JOB  CONTROL  LANGUAGE  FOR  EMPFIT 

The  following  JCL  is  necessary  to  execute  EMPFIT  on  the  HDL  IBM 
370/168  computer.  The  procedure  ANAFORT  is  used  to  supply  the  plotting 
software  while  the  object  version  of  EMPFIT  exists  on  permanent  file. 

//Job  Card 

/*JOBPARM  CARDS»20000 

//ST  EXEC  ANAFORT, PPELIB= ' HK3002 . EMPFIT' ,OUT=X,F3= 

//FORT. SYS IN  DD  * 

//LKED . SYS IN  DD  * 

INCLUDE  SYSLIB (EMPFIT) 

ENTRY  MAIN 

I//GO .FT10F001  DD  DISP= (NEW ,CATLG) ,VOL=SER=USER02 , 

//  UNIT=SYSDA,SPACE=(TRK, (2,2)) ,DSN=permanent  file 
name  for  curve  fit  output  data,DCB= (RECFM=VS) 

//GO.FTllFOOl  DD  DISP= (NEW.CATLG) ,VOL=SER=USER02 , 

//  UNIT=SYSDA,SPACE=(TRK, (2,2)) ,DSN=permanent  file 
name  for  derivative  output  data,DCB=(RECFM=VS) 

//GO.FT12F001  DD  DSN=permanent  file  name  of  disk 
file  input  data,DISP=SHR 
//GO. SYS IN  DD  * 


Input  Data 


// 

The  control  cards  in  braces  should  only  be  used  if  the  respective 
read  or  write  option  of  EMPFIT  is  employed. 


E-RAWAL  (X  101) 


BE6UESTE0  OPTIONS:  SCUBCE  .NCBBP.kOlBEF.CFT  (0 


•OPTIONS  IN  E F F E C T *SCUB  C E EBCDIC  NCLIST  NOCECB  OBJECT  NOPtP  NOFCBPBT  NOCOSTNT  NCXBEF  NOBIC  NOBNSF  NOT  EBP  FIBIIII 


LISTING  OF  EMPFIT  AND  SAMPLE  RUN  (cont'd) 
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SN  OOAA  AO  FORMAT I2E 10.31 

SN  00*S  SO  CUNTINUf 

SN  00*6  60  IFdCUT.EO.il  CALL  CUTOFF 

SN  00*8  RETURN 

SN  00*9  ENE 
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LISTING  OF  EMPFIT  AND  SAMPLE  RUN  (cont'd) 
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OPTIONS  IN  EFFECT:  NAME  I MA 1 N I NGOPTIMI2E  L INECCbNT  160  I SIZE (05  OOK ) AUTODBL (NONE  I 

SOUACE  EBCDIC  NOUST  NODECK  OBJECT  NONAP  NOFORMAT  NOGOSTMT  NOIRE  F NCAIC  N LAN  SF  NOTERM  FLAG)  II 
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LISTING  OF  EHPFIT  AND  SAMPLE  RUN  (cont'd) 
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Figure  B-l.  Example  of  data  preparation  by  using  disk  file  to  enter  data 
points,  using  time  cutoff  option,  deleting  data  points. 
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Figure  B-2.  Example  of  data  preparation  by  using  disk  file  to  enter  data 
points,  deleting  data  points,  adding  data  points. 
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IRBCXQINO  ?A0£  fti.ABy 


1 

3 3-1 .0001*03 

LEHP  PROBLEM  2 

42  4 1 . 000E-01  1.200E«08  5.000E«05  1 .000E*00  3.000E*09 

2 1 
1 .424 E-07-4.389E-05 
1 .563E-07-1  . 088  E-04 
1 .718E-07-2.948E-04 
1 .886E-07-6 .535 E-04 
2.072E-07-2.522E-03 
2 . 498 E- 07 -1.322 E -02 
2 . 745 E- 07 -2. 469 E -02 
3 .021 E-07-5.096E-02 
3.653E-07-4.974t-01 
4.003E-07-5.838E-01 
4 .380E-07-6.015E-01 
4 .793E-07-5.843E-01 
5 .259E-07-5.215E-01 
8 .368E-07-3.596E-01 
9.178E-07-3.372E-01 
1 .104E-06-3.185E-01 
1 . 326E- 06-3 .054  E -0 1 
1 .595  E-  06  -2  .884  E-0 1 
2.103E-06-2.789E-01 
2.529E-06-2.753E-01 

3 .041  £-06  -2  .64 1 E-0 1 

4 .399  £-06-1  .97  8 E-0 1 
5 .287E-06— 1 .459  E -0 1 
6.358E-06-9.011E-02 
8 .377E-06-2.672E-02 
1 .007E-05-8.379E-03 
1 .105E-05-4.926E-03 
1 .211E-05-3.658E-03 
1 .328  E-05— 3 . 954  E-0  3 
1 .4  56E-  05  -5 . 38  1 E-0  3 

1 .919E-05-1 .054E-02 

2 .307E-05-1  .324  E-02 
2 .773E-05-1  .464  E-02 
3.334E-05-1.533E-02 
3.656E-05-1.546E-02 
4 .009E-05-1 .531 E-02 
4.300E-05-1 .546E-02 
5.284E-05-1 .619E-02 
6 .353E-05-1 .480E-02 
7.638E-05-1.396E-02 
8 .0  80 £-05-1  .36  7E-02 
8 .374E-05-1  .33 7 E-02 
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Figure  B-3.  Example  of  data  preparation  by  using  cards  to  enter  data  points. 
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